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<N ■ Abstract 
> 

\o 

I Scattering problem by several bodies, small in comparison with the wavelength, 

' is reduced to linear algebraic systems of equations, in contrast to the usual reduction 

^ . to some integral equations. 

o 

'p^: 1 Introduction 

Acoustic or electromagnetic (EM) wave scattering by one or several bodies is usually 
studied by reducing the problem to solving some integral equations. In this paper we 
show that if the bodies are small in comparison with the wavelength, then the scattering 
problem can be reduced to solving linear algebraic systems with matrices whose elements 
^ ' have physical meaning. These elements are electrical capacitances or elements of electric 

■ and magnetic polarizability tensors. The author has derived analytical explicit formulas 

allowing one to calculate these quantities for bodies of arbitrary shapes with arbitrary 
desired accuracy (see P). 

We derive these linear algebraic systems and give formulas for the elements of the 
matrices of these systems. There is a large literature on wave scattering by small bodies, 
see and references therein. The theory was originated by Lord Rayleigh who 
understood that the main term in the scattered field is the dipole radiation if the body 
is small. Rayleigh did not give formulas for calculating the induced dipole moments 
for small bodies of arbitrary shapes. The dipole moments are uniquely defined by the 
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polarizability tensors. Therefore, the formulas, derived by the author (see [Q), allow one 
to calculate the dipole radiation for acoustic and EM wave scattering by small bodies of 
arbitrary shapes. 



2 Acoustic wave scattering by small bodies 

Let us start with acoustic wave scattering. Consider the problem 

m 

(A + e)u = in R3 \ ( U D^) (1) 

m=l 

u\s^=^, l<m<M, Sm.= dDm (2) 
u = Uq + v, (3) 

— ikv = o - , r := |a;| — > oo, (4j 

or \r J 

where A is the Laplacean, uq is an incident field which solves equation (P). Often, 
uq = e^^°''^, where a G 5"^ is a given vector and S"^ is the unit sphere. 
Let us look for the solution of the form 

u = UQ + y^ g{x,s)arr,{s)ds, g{x,y):=— (5) 

m=l 1^1, 

where 0"^, 1 < "nT- < M, are to be chosen so that the boundary conditions Q hold. The 
function (jSJ satisfies (P) and Q- (jH) for any am G L'^{Sm)- The scattering amplitude is: 



M 

Aia', a) = lim |x|e-*'=l^l = V — / 



e~''"''-'amds, a' := ^. (6) 



c \X\ 



Let 



and 



a := max diamDm, 

l<m<M 



d := mmdist{Dm, Dj] 



We assume 
Then 



A;a < 1, a < (i. (7) 
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so 



A{a\ a) = J2 4^ / ^n^ds := ^ e'^'^"'-^'", Q„ := / a„4s, (8) 

m=l '^'5™ m=l 

where Xm G Z^m and a' is defined in (0). Since is small, it does not matter which 
point Xm one takes in Dm- The Qm plays the role of the total charge on the surface Sm- 
If min^ |x — 3> a and Xm G D^, then 



u(a;) = Mo(2;) + r l + 0(fca + -) . (9) 

m=l ' ' 

Let us derive a formula for Q^. Using the boundary condition one gets: 



(10) 



where G ^m,. 

Since A;a -C 1, one has 



s) = go{sm, s) + 0{ka), 

where 

9o{s,t) : 



AtvIx — t\ 



Therefore equation is the equation for the electrostatic charge distribution am on 
the surface Sm of a perfect conductor Dm, charged to the potential 

Um ■= -Uo{Sm) " ^ g{Sm, Xj)Qj. 

The total charge on Sm is: 

Qm Cml^rrn 

where Cm is the electrical capacitance of the conductor with the shape Dm- The total 
charge is defined as: 

Qm - / (Tmds. 



Srf 



Therefore, one gets: 



Qm = Cmi -Uo{Sm) " ^ g{Sm, Xj)Qj j , 1 < m < M, (11) 

where Cm is the electrical capacitance of the perfect conductor with the boundary Sm- 
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Linear algebraic system (jllj) allows one to find Qj, 1 < j < M. If 

max V r < 1, (12) 

l<m<M ^ ATx\Sm — Xj\ 

then the matrix of the system has diagonally dominant elements and, consequently, 
can be solved by iterations. 

The approximate solution to the many-body scattering problem (^-(0} is given by 
formula ©, where Qm are determined from linear algebraic system (fTT|) . 

Let us give a formula from [1 for the capacitance of a perfect conductor D with the 
boundary S. Denote the area of S* by IS"!- We assume that the conductor is placed in the 
medium with the dielectric permittivity Sq = 1. In this case the approximate formula for 
the capacitance is (see [T], p. 26): 



C<»> = telSp N — / / — /..,/ ^»/'(«„-i.«„)*i (13) 




J Js Js "^st 

and the error estimate of formula ()13p is: 

-C\ = Oiq''), < g < 1, (14) 

where q depends on the geometry of S, and n = 1,2,3 is the approximation order. 

If the boundary condition 

un = Cu on Srn (15) 

is imposed in place of the Dirichlet condition (j2)), and C, is the impedance, then Cm in 
f|TT|) is replaced by 

C 

^'"^ i + cjc\s\)-'' ^^^^ 

see [T], p. 97. 
If 

UN |s„= 0, 1 < m < M, (17) 
then the formula for the solution to problem ((T)), (fTTj) . Q, (jH), is 



M 

U{x) = Uo{x) + ^ g(x, Xm)Vrn 
m=l 



A / \ _|_ \ ^ n .,du(^Xm) {x Xm), 



P 



p.,q=l 



where (x — Xm),v is the »-th coordinate of the vector x — x^, a ^ is the derivative with 
respect to the q-th. coordinate of x calculated at the point Xm, and /?pg,m is the magnetic 
polarizability tensor of Dm-, defined by the formula ([1], p. 98): 



S 



where Vm is the volume of Dm, the function a solves the equation 



a = Aa- 2Nq, 

N is the exterior unit normal to Sm, and 

f d 1 ^ ^ , 
= (^{t>dt, rst = \s-t\. 

The formulas for the tensor (3pq^m, analogous to the formulas (|T!?|) - (fT^ for the capac- 
itance, are derived in |Tl p. 55, formula (5.15)]. The unknown quantities Au{xm) and 
1 < m < M, 1 < g < 3, in p8|) can be found from the following linear algebraic 
system, analogous to (fTT|) : 

Au{xm) = Auo{xm) -k^ ^ g{xm,Xj)Vj[Au{xj) + ^ (3pg,jik^^^^^^ ^^"^ ^^''"'^ 



(19) 



du{Xm) _ dUo{Xm) ^ gff(Xn^,Xj) ^^ ^ .^^ du{Xj) {Xm - Xj),p . 

f^rf f^'y f't'y ^ J ^ ^^^^ J^Q^J I 'T* ^ I 

11 t-d j^m,j = l P,g=l JIM I J I 

(20) 

In (Hn)) we have used the equation 

^g{.x,y) = -k'^g{x,y), 

which holds if x y. 

From the linear algebraic system (fTU|) - (Pm) one finds the unknowns Au{xm) and 
1 < m < M, 1 < g < 3. 

If conditions ((Tj) hold, then system (fT^ -ipU j) has a unique solution which can be 
obtained by iterations. 

This completes the description of our method for solving many-body scattering prob- 
lem for small bodies and acoustic (scalar) waves. 
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3 Electromagnetic wave scattering by small bodies 



In the problem of electromagnetic (EM) wave scattering by many small bodies we assume 



a < A < d. 



[21] 



This assumption is more restrictive than ((7j). The reason is: in EM theory the fields 
are obtained by an application of first order differential operators, for instance Vx, to 
potentials, such as the vector potential. Applying this operator and calculating the field 
in the far zone one neglects the term | | compared with the term k. This means that 
the following inequality is assumed: 

1 1 



or 



d > A. 



In the acoustic wave theory the potential itself j^g{x^ s)ads has physical meaning, it is 
the acoustic pressure, and this pressure is studied. Therefore, the condition X does 
not appear. 

Condition ((7j) allows one to have many small particles on the distance of order A, 
while condition PHI . namely the inequality d ^ X, does not allow this. Recall that d is 
the minimal distance between two neighboring particles. The formula for the scattering 
amplitude, analogous to (jH)), for EM wave scattering by small bodies is (see [2]): 



1 



(22) 



m=l 



Here U = {^) is a 6-component vector, Sm is a 6x6 matrix, the scattering matrix, Eq and 
/io are dielectrical and magnetic parameters of the medium, in which the body Dm is 
placed, and 9, 9' are the unit vectors in the direction of the incident and scattered waves, 
respectively. These vectors were denoted a and a' in Section |21 We have changed the 
notations because in EM theory a denotes the polarizability tensor. 
The formula for S is (cf. [2]) 



Sr, 



^aE - 9'{9', aE) 



\ 



MP-0'i9',PH))j 



(23) 



Here Vm is the volume of Dm, 



a is the electric polarizability tensor of D^, P is the 
magnetic polarizability tensor of D^.- In pp. 54-55] the author derives analytical 
formulas for calculation of the polarizability tensors a and (3, 



/i + ^0' 



7 ■-- 



e + Eo 



(24) 
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Tensors (3 and (3 are expressed through the polarizabihty tensor a = 0(7). One has 
7 = —1 if £ = 0. Here [■, ■] is the vector product, (■, ■) is the scalar product. 

The analytic formula from p. 54, formula (5.9), for the tensor a = aij{'y), 1 < i, 
j < 3, that we referred to above, is analogous to formulas (fT^ - (fT^ for the electrical 
capacitance. The incident direction 9 enters via the vectors E and H, which depend 
on 9. These vectors are calculated in formula (j!23|) at the point Xm- The values of these 
vectors are determined from a linear algebraic system of equations. This system is derived 
similarly to the derivation of the systems (fTT|) and (jiyp - ipUj) . We do not write down this 
system since it would take much space, but the ideas are the same as the ones used in 
the derivations of (dH) and (fT!I|l -CT. 

4 Conclusions 

In this paper it is shown how to reduce rigorously the many-body scattering problem 
to linear algebraic system in the case when the bodies are small in comparison with the 
wavelength. The theory is constructed for acoustic and EM wave scattering. The basic 
physical assumptions are ((Tj) for acoustic scattering, and (PT|) for EM scattering. 
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